Confinement Models at Finite Temperature and Density 
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In-medium chiral symmetry breaking in confining potential models ol QCD is examined. Past 
attempts to analyse these models have been hampered by infrared divergences that appear at non- 
zero temperature. We argue that previous attempts to circumvent this problem are not satisfactory 
and demonstrate a simple resolution. We also show that the expectation that confining models do not 
exhibit a chiral phase transition is incorrect. The effect of summing ring diagrams is investigated and 
we present the first determination of the temperature-density phase diagram for three model systems. 
We find that observables and the phase structure of the confinement models depend strongly on 
whether vacuum polarisation is accounted for. Finally, it appears that standard confinement models 
cannot adequately describe both hadron phenomenology and in-medium properties of QCD. 
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I. INTRODUCTION 

The properties of QCD at finite temperature and den- 
sity find applications in topics as diverse as the nature 
of proto-neutron stars, early universe cosmology, and ex- 
periment at RHIC and the LHC[l[. Unfortunately, many 
of the properties of interest are nonperturbative, and, 
with the exception of lattice techniques, tools for dealing 
with nonperturbative field theory remain rudimentary. 
Heightening our discomfort is the continuing statistical 
'minus sign' problem in lattice field theory that is present 
at finite chemical potential [H. Furthermore, well-known 
problems of infrared (IR) divergences at high order in 
perturbation theory persist Q. Even old hopes that QCD 
at large temperature is perturbative may be misplaced 
since large temperature QCD corresponds to the dimen- 
sionally reduced QCD of three dimensions, which is also 
confining. Thus it is of interest to develop nonperturba- 
tive methods for analysing QCD and model systems in 
medium. 

Here we examine the properties of two simple models 
of confinement that are motivated by QCD in Coulomb 
gauge. The analysis of such models dates from the mid- 
1980s and parallels extensive work with Nambu-Jona- 
Lasinio (NJL) models (We note that NJL models are 
not renormalisable and do not exhibit confinement.) In 
the following we will examine a contact model in detail, 
focussing on coupling constant dependence of the phase 
diagram, possible critical points, and the effects of incor- 
porating polarisation in the formalism. We demonstrate 
that a surprisingly rich structure emerges. 

We also examine confining potential models with the 
intent of clarifying several points in the literature. The 
first of these concerns the existence of infrared diver- 
gences in the temperature-dependent gap equations that 
naively renders them useless. The second issue is whether 
a linear model can support a chiral phase transition. We 
resolve these issues and also examine the effects of po- 
larisation on the phase diagram. As far as we are aware 
this investigation presents the first determination of the 
phase structure of confining models and the first exami- 



nation of general vacuum polarisation effects in contact 
and confining models. 



II. CONFINEMENT MODELS AND INFRARED 
DIVERGENCES 

A. Linear and Contact Models 

Our starting point is the Hamiltonian of QCD in 
Coulomb gauge. Coulomb gauge is especially appropriate 
for the study of in-medium properties of QCD because all 
its degrees of freedom are physical. The usual demerit as- 
sociated with non-manifest covariance is obviated by the 
presence of the heat bath. Upon neglecting transverse 
gluons, the QCD Hamiltonian takes the form 



^(-i7-V+m)f + - p a (x)V(x-y)p a (y) (1) 



H 



where p a = ip' ! T a ip is the colour quark current and T a is 
a generator of SU(N). 

Neglecting transverse gluons ruins the gauge and 
Lorentz invariance properties to the theory. However, it 
has been argued that (dynamically) massive transverse 
gluons provide a more accurate basis for the exploration 
of low energy properties of QCD[j| Q. Thus, as long 
as the induced error is not too great, gauge invariance 
should be approximately respected. Similarly, the full 
theory is Lorentz invariant, even if it is not manifestly 
covariant. Accurate truncations should accurately pre- 
serve this property. 

The instantaneous interaction kernel corresponds to 
the Coulomb potential in QED; in the case of QCD it 
can be written as the vacuum expectation value of the 
Coulomb operator 
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Here D ab = VS ab + gf acb A c is the adjoint covariant 
derivative, a, b are colour indices, and f2 is the full vac- 
uum. This interaction is often modelled as a linear con- 
finement potential: 



V(r) 



-br. 



V(q) 
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The string tension is denoted b and its phenomenological 
value is approximately 0.2 GeV 2 . 

An alternative that matches to perturbation theory is 
provided by the Richardson potential, 



V(q) 



-in 



4 gWog(l + g 2 /A 2 ) 
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with (3q = 11 — A 2 — 26/3o, and n/ is the number of 
quark flavours. 

Finally, we shall consider a simple contact model de- 
fined by 



employ the Schwinger-Dyson equations as our organising 
principle. 

The Schwinger-Dyson equation for the full fermion 
propagator in potential models such as those employed 
here is shown in Fig. [1] The dashed line in this figure rep- 
resents an application of the instantaneous interaction. 
Notice that this line is not dressed and does not have an 
analogous Schwinger-Dyson equation because it is not 
dynamical. Inserting the tree level approximation to the 
fermion four-point vertex of Fig. Q] yields the equation 
represented in Fig. [2] This equation represents a series 
of diagrams that yield the rainbow-ladder approximation 
to the fermion propagator, corrections to the rainbow- 
ladder approximation generated by a dressed vacuum po- 
larisation insertion, and corrections due to vertex dress- 
ing. 
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V (r) = —8(f). 



(5) 



The scale A is introduced to set dimensions and will be 
used as an ultraviolet cutoff in this model. Of course 
the contact model is not confining, however it permits 
detailed analysis of in-medium effects since the resulting 
equations are considerably simplified. It also serves as a 
reference point for the confining models. 
The partition function is defined as 



Z[r),v\ 



Dip Dip exp[— A + rjip + iprj\ 



(0) 



with 



A = I dTd 3 xip(~f (d T - n) -ij -V + m)ip + 



dTd i xdT , d i yp a {x)V(x-y)5(T-T l )p a (y) (7) 



Notice that a quark chemical potential term, propor- 
tional to fi, has been added to the action. 

We employ the imaginary time formalism for evaluat- 
ing the partition function. In particular, the time inte- 
gral reduces to sums over bosonic and fermionic Matsub- 
ara frequencies. These sums are evaluated by integrating 
over appropriate contours using the Orsay representation 
of the propagators. 



B. Schwinger-Dyson Equations 

Our goal is to study the in-medium properties of chi- 
ral symmetry breaking in the models defined in section 
III Al Thus nonperturbative methods are required and we 



FIG. 1: Schwinger-Dyson Equation for the full fermion propa- 
gator in potential models. Minus signs are not made explicit. 
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FIG. 2: Expanded Schwinger-Dyson Equation for the full 
fermion propagator. 

If one neglects vertex correction diagrams, it is possi- 
ble to sum all dressed vacuum polarisation insertions by 
rewriting the equation represented in Fig. [5] as shown in 
Fig. [31 These coupled equations, called the gap equa- 
tions, form the starting point for our investigation of 
dynamical mass generation. It is important that the 
vacuum polarisation fermion loop of the second equa- 
tion utilises dressed fermion propagators. Failure to do 
so would yield a branch cut that signifies the decay of 
dressed fermions to bare fermions - which is clearly phys- 
ically unreasonable. 
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FIG. 3: Gap Equations: summing polarisation insertions in 
the truncated Schwinger-Dyson equations. 

The second gap equation (Fig. [3]) implements the "ring 
approximation" to the full interaction potential. This 
approximation was first employed by Brueckner and Gcll- 
MannQ to solve the infrared divergence problem in the 
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electron self energy of the degenerate electron gas. As we 
shall discuss shortly, it will serve the same purpose here. 
We define the polarisation as 



d 3 P 

(2tt) 3 



ti[ 10 S(khoS(p + k)} (8) 



where n/ is the number of quarks (considered of equal 
mass in the following), S is the full quark propagator 
(discussed more fully below), and k is a four vector de- 
fined as 



k^ = (iu n + fJ,,k), 



(9) 



the fermionic Matsubara frequency is given by w n = 
(2n + 1)tt//3, and \i is the quark chemical potential. The 
colour trace yields the factor of \ in Eq. [H The expres- 
sion for the polarisation is usually split into vacuum and 
matter components, with the vacuum contribution being 
defined as H(k:T — > 0, /i — > 0). Rcnormalisation follows 
the standard vacuum formalism and hence only affects 



n v 



With this definition the ring potential is given by 



Vrin g {qo,q) = 



V(q) 



l-IL(q ,q)VW 



(10) 



To establish contact with well-known results we note 
that 



lim U(p Q = 0,p) 

p-s-0 
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in the case of a massless bare quark. Of course the full 
computation must be made with dressed quark propaga- 
tors and hence forms part of the coupled gap equations. 

Incorporating polarisation effects in the gap equations 
can be very important. For example, quantum electrody- 
namics in three dimensions is a (logarithmically) confin- 
ing theory; however, including polarisation effects due to 
massless fermions completely screens the confinement po- 
tential, leaving a Coulombic heavy fermion interaction [9]. 
Similarly, increasing the number of quarks in QCD even- 
tually drives the theory into a conformal window with no 
confinement [l(| ■ 



The scalars fi, A, and B are functions of ko and |A;|. 
Note that VJ-i ng depends on ko and k, substantially com- 
plicating the solution to the gap equations. However, 
the dominant infrared contribution to the ring potential 
is obtained in the static limit, ko — > and we employ 
this limit in the following. Under these conditions the 
gap equations become: 



A(p) 



2 J (2tt) 



d 3 q T/ ,„ Aqp-q 
3 King (p - q) — — @{q) 



E q p 2 
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q)[n{q) - n{q)} 



(13) 



We have introduced the colour factor Cf = (N 2 — 
1)/(27V). The thermodynamic function is defined in 
terms of the quark densities as 



0(q) = 1 - n(q) - n{q) 



with 



and 
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n{p) 
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exp(/3(^ p +/2)) + 1' 



(14) 



(15) 
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Similar equations have been considered by Kocic [ll| 
and more recently in Refs. [H, EH- Kocic presents an- 
alytic results for a contact model very similar to ours, 
but without consideration of vacuum polarisation effects. 
All three papers consider the linear confinement model 
at zero temperature; Ref. [l2| incorporates ring correc- 
tions and examines the effect on charmonium dissocia- 
tion, but neglects the vacuum part of the polarisation 
without comment. A contact model is considered at fi- 
nite temperature and density in Ref. 03; additional 
related studies are listed in Ref. 



In the case of the contact potential or the ring con- 
tact potential in the static and long wavelength limits 
(discussed below) these equations simplify further to 



C. The Gap Equations 

The diagrams of Fig. [3] represent four coupled integral 
equations; three involve the scalar functions defined in 
the general expression for the in-medium inverse quark 
propagator: 



S 1 (fc) = i(uj n — — 7 • kA — B. 



A 
B 

Elip) 



1, 
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Cf 
2 J (2tt^ 

C F f d 3 q 
2 J (2tt) 
S 2 , 



d3q ^(o)lteo(g), 



3 VS£@)[n(q)-n(q)] 



(12) 
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(18) 
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where Qq is defined in terms of -Bo- 
More generally, the ring contact case is as complicated 
as the linear model since the interaction becomes 



t rcon 
*ring 



(P) = 



A 2 i 



\ n(o,p) 



(19) 



Since the chief role of the polarisation function is to regu- 
late the infrared divergence that appears in the gap equa- 
tions (discussed in the next section), it is appropriate to 
consider the approximation: 



Vring(g) 



V(q) 



l-U(0,q^0)V(q) 



(20) 



Indeed, as Eq. [501 illustrates, in the linear case V(q) 
dominates Vri ng for q 2 3> m g yb and it is reasonable to use 
the zero momentum limit of the polarisation function in 
the gap equations. We shall confirm this in the following. 
In this limit the vacuum contribution to the polarisation 
function vanishes and we need only consider the matter 
contribution. In fact no ultraviolet divergences remain in 
the linear model and we do not consider renormalisation. 
Of course the contact model must be cutoff at the scale 
A, as previously indicated. 

In the contact case, this approximation implies that 
Vri n g is a constant and hence B is a constant. Thus the 
integral equation for B in the contact case becomes a 
relatively simple algebraic equation, greatly simplifying 
the problem. 

In the low momentum limit the explicit expression for 
II becomes 



n(0,g^0) = -4n / /3 / ^[2 + + e"^]- 1 

J {^) 

(21) 

where E(k) depends on A and B and is naively infrared 
divergent (cf. Eqs. Upland section u*TD|) . 

Finally, the dynamical quark mass is given by 



M(p) 



B(P) 
A(p) 



(22) 



and a gap equation can be derived for this function di- 
rectly: 



M{p) — m H — — 



Cp f d 3 g 
M(q) M(p)p-q 



E (q) E (q) p 2 



V riDg (p- q) 



e( 9 ) 



(23) 



where now Eo(q) = \/q 2 + M(q) 2 . It is important, how- 
ever, to note that still depends on the full single par- 
ticle energy, E(q). 



D. Infrared Behaviour of the Gap Equations 

In the absence of polarisation effects the linear model 
gap equations of Eq. [T3] are naively infrared divergent. 
Specifically as q — > p + S on the right hand side, the 
potential diverges as V(S) ~ 8~ A . Thus the quark en- 
ergy E(q) is divergent, the thermodynamic function ap- 
proaches unity, and temperature effects disappear from 
the in-medium gap equations. We stress that it is not 
possible to write the equations in terms of IR-finite quan- 
tities, in contrast to the zero temperature case. This non- 
sensical result was first noted by Davis and Matheson[l6| 
who suggested that one should make the ad hoc replace- 
ment 



E(q) -+ E(q) - E(0) 



(24) 



in the expression for 0(g). This removes the infrared 
divergence and yields a sensible in-medium gap equation. 

An alternative suggestion was made by Alkofer et 
«L[17I]. who simply replaced E(q) with Eo(q). The utility 
of this Ansatz is that the equation for M is infrared finite 
as evidenced in Eq. [531 We shall denote this procedure 
as 'AAL' in the following. 

Finally, the Orsay group has argued that infinite quark 
energies are physically reasonable and reflect the absence 
of individual quarks in the physical spectrum [18j. They 
therefore reformulate the thermodynamic trace to sum 
over colour singlet states only. This greatly complicates 
the thermodynamic trace but eventually yields a simple 
result: momentum space integrals such as appearing in 
Eqs. [T3] are restricted to q ^ p, thereby eliminating the 
infrared divergence. 

None of these resolutions seem appropriate. The sub- 
stitutions of Refs. [HI and [l?} are completely ad hoc. 
However, it is possible that they are reasonable approx- 
imations, and we investigate this in the next section. 
The Orsay approach is physically motivated, however, it 
should not be necessary to explicitly remove nonsinglet 
degrees of freedom from the thermodynamic trace be- 
cause they are automatically removed by the Boltzmann 
factor. 

We note the following: (i) infrared divergences are 
normally removed by considering additional physical ef- 
fects (such as bremsstrahlung) or by summing additional 
classes of diagrams (as with the resolution of the IR di- 
vergence problem in the degenerate electron gas[8[ or in 
the hard thermal loop (HTL) partial resolution ofgauge- 
variance and other problems in thermal QCD[19|]) (ii) 
'artificial' infrared divergences sometimes appear in prob- 
lems involving confinement potentials 20] . 

We examine both issues in the following. Specifically, 
issue (i) has been addressed by summing the ring dia- 
grams in the gap equations. This has the effect of remov- 
ing the infrared divergence in all expressions, as evident 
in Eq. 1201 Surprisingly, issue (ii) is also relevant here. 
Indeed, the gap equations are trivially IR regulated by 
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subtracting zero. For example, the third of Eqs. [13] can 
be rewritten as 



Bm 



m - 



d 3 q 

(27T) 



E. 



(25) 

While this yields a formally IR-finite equation, we find 
that it is still strongly dependent on the infrared regime, 
thus it remains useful to employ the ring potential in the 
formalism when performing numerical computations. 

We examine the regulated gap equations with the bare 
and ring potentials in the following section. 



III. PHASE STRUCTURE 

Our results will be presented as plots of the dynamical 
mass at zero momentum as a function of chemical po- 
tential and temperature. Thus we determine the phase 
structure of the contact and confinement models. As 
far as we know these are the first computations of these 
phase diagrams for confining potentials. The results for 
all potentials with ring corrections are also new. 

The coupled gap equations are solved with a variety 
of techniques. Our typical approach is to compute on a 
momentum grid. Unknown functions A, B, and 14ing are 
determined by minimizing an appropriate functional that 
represents the gap equations using a modified Levenburg- 
Marquardt algorithm. Alternatively, we attempt itera- 
tive solutions (which are often not stable) and combina- 
tions of iteration and minimisation. We also employ a 
variety of analytic approximations to quantities such as 
M and V r i ng and confirm the more complete results. For 
example, computing the full momentum dependence of 
II is time-consuming and it is useful to expand it as a 
Taylor series in momentum. The results presented here 
represent several months of single core computer time. 



A. Contact Model Results 

The contact model is defined in terms of the coupling A 
and ultraviolet cutoff A. It is IR-finite and hence the only 
issues here are the effect of the ring interaction and, of 
course, the new phase structure shown below. Contact 
models such as ours have a critical coupling A c , below 
which no chiral symmetry breaking occurs. Since this 
is a zero-temperature property of the model and n mat 
vanishes in this limit, the critical coupling is the same 
for the bare and ring models: 



A r = 6?H 



(26) 



To set conventions we choose to work in the bro- 
ken phase, setting A = 1.5A C and then fixing the scale 
A = 370 MeV to yield reasonable approximations for the 



dynamical quark mass and chiral restoration tempera- 
ture, M{T = 0, fi = 0) « 260 MeV and T c « 150 MeV. 
The chiral condensate is given by the expression 



(MO = —j 



k Z dk 



B(k) 
E(k) 



(27) 



and is numerically (-150 MeV) 3 , which is approximately 
a factor of two too small (in linear dimension). Finally 
the critical chemical potential is approximately 300 MeV. 
The numerical values of these observables shift slightly 
with the coupling, for example for A = 3A C one finds 
A = 140 MeV, T c = 135 MeV, and ^ c = 250 MeV. 

It is the fact that the zero temperature gap equation 
has a critical coupling that drives the finite temperature 
and density phase transitions. Specifically, the presence 
of the thermodynamic function O at T > can be mod- 
elled as an effective coupling: 



A(T, M ) = A<e(?;T )M )) 



(28) 



where the angle brackets denote some average over the 
integrand of the gap equation. Once A(T, /i) drops be- 
low A c the model makes the transition to the symmetric 
phase. 

The phase diagram for the bare contact model is shown 
in Fig. SI In the following we will refer to Mq defined 
as M(k — > 0;T, /i). This represents a 'constituent quark 
mass' and is a useful and simple proxy for the chiral con- 
densate. We note that the simplicity of the contact gap 
equations permit a detailed analysis of the phase struc- 
ture, including a precise determination of the order of 
the phase transition. This is useful since, as the figure 
illustrates, the phase structure is quite complex. We find 
that the chiral symmetry restoration phase transition is 
second order for all chemical potentials below a critical 
density 



^*(A = 1.5A C ) w 0.53A. 



(29) 



For higher chemical potential the transition is first or- 
der. Furthermore, a second solution to the gap equations 
develops - these are the lower lines seen in the figure 
between /x* and = 0.9. These solutions indicate the 
presence of a first order phase transition; however, they 
are unphysical in that they have higher free energy, as 
demonstrated by their lack of continuity with the low-/i 
solutions. 

Including the zero-momentum vacuum polarisation 
function (with n/ = 1) in the gap equations induces non- 
trivial dynamical mass, temperature, and chemical po- 
tential dependence in the kernel of the gap equation. The 
results are presented in Fig. [5j It is apparent that the ad- 
ditional mass dependence causes the bifurcated solution 
to exist for all chemical potential. The phase transition is 
first order for all values of temperature and chemical po- 
tential. Furthermore, the numerical value of the critical 
temperature is strongly affected: 
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FIG. 4: Dynamical mass vs. temperature and density for the 
contact model (A = 1.5A C , all quantities in units of A). 



T c (fi = 0; A = 1.5A C ; bare) = 0.38A (30) 

drops to 

T c (n = 0; A = 1.5A C ; ring) = 0.17A. (31) 

For A = 3 A c the analogous results are T c (bare) = 0.95 A 
and Turing) = 0.30A. Of course, ^ c does not change 
because it is defined at T = and II = at this point. 
Thus the shape of the critical region and the phase tran- 
sition are strongly affected by moving beyond the bare 
interaction approximation in the gap equation. 




FIG. 5: Dynamical mass vs. temperature and density for the 
contact model with the static, zero momentum ring potential 
(A = 1.5A C , all quantities in units of A). 



Finally, we consider the case where the full momen- 
tum dependence is retained in the vacuum polarisation, 
leading to a momentum dependent ring potential. The 
results are shown in Fig. [5] We see that the dynamical 
fermion mass has dropped to M = 0.23A, as have the 
critical temperature and chemical potential T c = 0.06A, 
and fj, c = 0.22A. However, upon rescaling A we find 
that all of these quantities remain essentially unchanged. 
Furthermore, it appears that the phase transition is first 
order throughout the phase diagram, as it is in the zero- 
momentum ring case (determining the order of the phase 
transition is difficult with momentum dependent inter- 
actions because algebraic methods are not applicable). 
Thus it appears that including the full momentum de- 
pendence of the polarisation function and ring potential 
alters the numerical values of quantities in units of the 
cutoff, but leaves physical quantities approximately in- 
variant. 




FIG. 6: Dynamical mass vs. temperature and density for the 
contact model with the static ring potential (n mat + II vac , 
A = 1.5A C , all quantities in units of A). 

We have also studied the A-dependence of the phase 
diagram for the 'bare' contact model. We find that all 
quantities scale roughly linearly with A. For example, to 
a good approximation 

M (A) = 0.011 ■ AA- (l - e -o.o28(A-A c )^ ( 32 ) 

Furthermore, the relative strengths of the dimensional 
quantities seem to be roughly independent of the cou- 
pling, and we find: 

Mc(A) ps M (A) « 2T C (A) « 2/x* (A) (33) 

The coupling dependence of the phase diagram for the 
bare interaction case is illustrated in Fig. [7J The approx- 
imate linear scaling is evident. This is also true for the 
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critical point /i* down to A ss 1.2A C . Below this point 
the hrst order phase transition disappears and the entire 
diagram represents a second order phase transition. 
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FIG. 7: The bare contact interaction phase diagram for a 
variety of couplings (£ = A/A c ). 



B. Confinement Model Results 

The confinement models exhibit the full IR singulari- 
ties that have driven the discussion. Unfortunately, the 
models are much more difficult to analyse and one must 
rely on accurate numerics to establish the order of phase 
transitions. The relation of the chiral restoration and de- 
confinement transitions is also of interest. Longstanding 
lattice gauge theory results indicate that these transition 
temperatures coincide [2lj|. Since the confinement poten- 
tial is fixed bare models cannot reproduce this behaviour. 
Indeed Davis and Matheson have argued[lj| that all con- 
fining potentials break chiral symmetry and since the 
only in-medium effect is the addition of the factor in 
the gap equation kernels, it is feasible that there can be 
no transition to the chirally symmetric phase. Never- 
theless, Alkofer et al. have numerically found a phase 
transition. 

As shown in Fig. [5] we confirm the existence of a phase 
transition in the confinement model with the AAL pre- 
scription (and we extend their results to the T—\i plane). 
We speculate that the argument of Davis and Mathe- 
son has failed because for high temperature and low mo- 
mentum the thermodynamic function is approximately 
<d(q,T,fi — 0) ~ q/T, and hence the effective confining 
potential is replaced by V ~ q~ 4 —> q~ 3 , which is not 
confining. Perhaps this is sufficient to drive the observed 
behaviour. 

Numerical values for the linear potential are T c 38 
MeV, fi c « 75 MeV, M rj 80 MeV, and the conden- 
sate is determined to be (— (i/j?/;)) 1 / 3 ps 110 MeV. Note 
that we interpret the dynamical mass as a constituent 
quark mass. All of these quantities fall below their phe- 
nomenologically expected values; however, using a string 



tension of about 1.8 GeV 2 brings them all into reason- 
able agreement. In fact, Alkofer et al. employed such a 
large string tension in their computation. This is unfortu- 
nately at odds with well established hadron phenomenol- 
ogy; a point which we will discuss further below. 

The chiral restoration transition line is plotted in Fig. 
[5] A clear inflection point is visible at (T+, sa (20, 43) 
MeV. We identify the right side of this point as a region 
of first order phase transitions, while the left side is sec- 
ond order. This identification is based on the continuity 
of M displayed on Fig. [31 the evident first order phase 
transition seen at T = 0, the relative lack of stability of 
the solution algorithm for large /i, and the appearance of 
secondary solutions (not shown) for // > /i*. The exis- 
tence of a tricritical point is in keeping with expectations 
forQCD@. 
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FIG. 8: Dynamical mass vs. temperature and density for the 
bare linear AAL confinement model. 
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FIG. 9: The bare linear AAL phase diagram. 

The results of Fig. [8] have been computed with the 
ad hoc AAL prescription. We now consider modifying 
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the infrared behaviour of the quark interaction with the 
vacuum polarisation diagram of Fig. [3] Our first study 
is in the static, low momentum limit with n/ = 1. In this 
case the vacuum contribution to II vanishes and we need 
not consider renormalisation of the model. 

Our results are shown in Fig. [TO] The AAL and ring 
results must agree at the origin. The surfaces at /i = 
look very similar, however the critical temperature has 
dropped from 38 MeV to 10 MeV. The critical chemi- 
cal potential remains fixed at /z c ~ 75 MeV. Thus the 
shape of the phase diagram has been severely distorted 
by incorporating vacuum polarisation in the model. Fur- 
thermore, as with the contact model, it appears that all 
phase transitions are now first order, and the tricritical 
point no longer exists. Thus, as anticipated, including 
dynamical quarks in the model can have a dramatic ef- 
fect on thermodynamic properties. 




FIG. 10: Dynamical mass vs. temperature and density for 
the linear static long wavelength ring approximation 

The vacuum polarisation function introduces explicit 
temperature and density dependence to the quark in- 
teraction, which raises the possibility of explicit quark 
deconfinement in the model. It is possible that this de- 
pendence causes the potential to deconfine at a critical 
temperature. However, it is more likely that the potential 
deconfines for all nonzero temperature. Indeed, in per- 
turbation theory one can approximate the ring potential 
as 



V xins (q,T^ = 0) 



Qitb 



irbT 2 ' 



(34) 



The Fourier transform of this potential is linear when 
T = 0. When T > the potential is linear at small 
distances, has a transition region at r ~ (ttW 12 )- 1 / 4 , and 
approaches zero at large distances, so that deconfinement 
is natural, although not sudden. A careful analysis of 
deconfinement awaits a study of QCD. 



Finally, we consider incorporating the full momentum 
dependence of the polarisation function in the gap equa- 
tions. It is important to note that the vacuum contri- 
bution to the polarisation function is not zero for non- 
zero momentum. This contribution must be renormalised 
at the expense of introducing another parameter to the 
model. In practice, the relatively good agreement with 
heavy quark spectroscopy would be ruined by this pro- 
cedure. We therefore consider II vac to already be incor- 
porated in the model interaction, and simply consider 
matter contributions to the ring potential. 

Unfortunately, the numerical solution of the gap equa- 
tions is substantially complicated in this case. Solving 
the full set of coupled integral equations can be time con- 
suming and unstable. We have found that setting A(q) 
to be large provides a very good approximation to the 
full gap equations, and this approximation has been em- 
ployed for most of the linear ring results presented here. 
An additional complication is provided by the require- 
ment to compute the polarisation at many grid points at 
each step of the computation. However, we have found 
that expanding H(q) to 0(q 4 ) is very accurate, and this 
speeds the solution tremendously. 

In the end, the results are very similar to those of the 
zero momentum case, as anticipated above. For example, 
the dynamical mass M(k — 0,T = 0,[i — 0) is lowered 
by approximately 1 MeV when incorporating the momen- 
tum dependence of the vacuum polarisation. Similarly, 
the critical temperature remains the same to the accu- 
racy we compute. 

The results presented in this section are for the linear 
potential of Eq. [3J We have performed the same compu- 
tations with the Richardson confinement model of Eq. 2] 
This permits examining the dependence of our results on 
the details of the confinement model. In particular, the 
Richardson potential incorporates a running Coulomb in- 
teraction, and it has been speculated that this can en- 
hance the chiral condensate. We find, however, that the 
Richardson potential yields nearly identical results as the 
linear potential. For example, the AAL Richardson dy- 
namical mass is only 1 MeV less than the corresponding 
linear mass. Similarly, the critical temperature drops 0.7 
MeV and the critical density raises 3.2 MeV. 



IV. CONCLUSIONS 

We have examined bare and ring versions of a contact 
model that mimics the structure of QCD in Coulomb 
gauge. First results for the phase diagram of the ring 
contact model are presented here. The bare model ex- 
hibits interesting behaviour, including approximate lin- 
ear scaling of T c , fi c , and the dynamical mass with the 
coupling. The numerical values of these quantities can be 
placed in rough agreement with QCD expectations with 
an appropriate choice of the cutoff scale. Surprisingly, 
we find a tricritical point for all models with A > 1.2A C , 
again in keeping with expectations for QCD. However, in- 
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corporating ring-type Schwinger-Dyson equations in the 
formalism changes these conclusions dramatically. Ring 
values for /i c , Mq and the condensate do not change sub- 
stantially, however, the critical temperature drops by a 
factor of three, ruining the phenomenology. Further- 
more, the region of second order phase transitions dis- 
appears and only a first order phase transition appears 
in the phase diagram. Given the strong n/-dependence 
expected in gauge theories, perhaps the strong effects 
seen here should not be surprising. 

For the bare linear and Richardson models with the 
AAL infrared prescription we confirm the existence of a 
second order phase transition at small chemical poten- 
tial. For chemical potential larger than /z* s» 43 MeV the 
phase transition becomes first order. The appearance 
of any phase transition is somewhat surprising, since it 
is in conflict with the reasonable expectations of Davis 
and Matheson. The numerical values for the dynami- 
cal mass, chiral restoration temperature and density, and 
chiral condensate are all in agreement with QCD expec- 
tations if the string tension is increased to a value of 
1.8 GeV 2 . Unfortunately, this is in severe conflict with 
well-established quark model phenomenology and lattice 
gauge results that require a string tension of approxi- 
mately 0.2 GeV 2 . It is thus apparent that the simplest 
confinement models cannot both reproduce thermody- 
namic and spectroscopic quantities with any reliability. 
Of course this conclusion depends on the approximations 
we have made. However, the large discrepancy seems 
difficult to overcome and we expect that simple confine- 
ment models are incapable of describing in-medium prop- 
erties of QCD. These conclusions hold for both linear and 
Richardson models. 

We have noted that two resolutions to the confinement 
model infrared problem exist: a direct physical resolution 
is to employ the ring approximation and a simple math- 
ematical resolution involves the proper regularisation of 
the Fourier transform of the linear potential. We find 
sufficient numerical infrared instability to warrant em- 
ploying the physically reasonable, and presumably more 
accurate, ring Schwinger-Dyson equations. Results in the 
static/low momentum and static/finite momentum lim- 
its were found to coincide very well. As with the contact 



interaction, ring diagrams induce strong effects in ther- 
modynamic observables, causing the critical temperature 
to drop by a factor of four and changing all phase tran- 
sitions to first order. 

Pisarski and McLerran have recently used the large 
N c limit to argue that a 'quarkyonic' QCD phase may 
exist [22| . This phase is proposed to exist at high chemical 
potential and low temperature and to be characterised by 
chirally symmetric but confined quark matter. (We note 
that the argument is supported by observations that chi- 
ral symmetry breaking should be irrelevant high in the 
zero temperature hadron spectrum [23| . ) The large N c 
limit suppresses quark loop effects, thereby yielding an 
interaction that is independent of the chemical potential. 
Thus the linear no-ring model presented here can be con- 
sidered an implementation of the large N c scenario. Our 
results then confirm the idea that a confining but chi- 
rally symmetric phase can exist. Of course, decreasing 
the number of colours implies an increasing importance 
of vacuum polarisation with commensurate changes in 
phenomenology, as noted above. A full assessment of the 
quarkyonic matter scenario must await a nonperturba- 
tive investigation of more realistic models of QCD. 

The application of finite temperature and density 
Schwinger-Dyson methods to the QCD Hamiltonian will 
be of great interest. As illustrated here, one of the over- 
riding considerations will be the construction of a ro- 
bust truncation scheme for Hamiltonians with infrared 
enhanced interactions. Indeed, the Ward identities relate 
vertices and propagators, so one expects that if correc- 
tions to propagators are necessary, one should truncate 
full vertices very carefully. Of course this observation is 
borne out by the experience with the effective field theory 
generating the HTL limit of QCD. 
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